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The similarity of the mathematical description of random-field spin systems to orthogonal 
frequency-division multiplexing (OFDM) scheme for wireless communication is exploited in an 
intercarrier-interference (ICI) canceller used in the demodulation of OFDM. The translational sym- 
metry in the Fourier domain generically concentrates the major contribution of ICI from each sub- 
carrier in the subcarrier's neighborhood. This observation in conjunction with mean field approach 
leads to a development of an ICI canceller whose necessary cost of computation scales linearly with 
respect to the number of subcarriers. It is also shown that the dynamics of the mean-field canceller 
are well captured by a discrete map of a single macroscopic variable, without taking the spatial and 
time correlations of estimated variables into account. 

PACS numbers: 05.40.-a, 75.10.Hk, 84.40. Ua, 88.80.ht 



I. INTRODUCTION 

Wireless communication technologies play a significant 
role in the modern information society. As of the end 
of 2010, there are more than 4.6 billion mobile-cellular 
subscriptions in the world [l[ , and the use of wireless de- 
vices (such as personal digital assistants and GPS units) 
is ever increasing. To keep up with the accompanying 
rapid growth in data traffic, communication efficiency 
of today's wireless communication systems must thus be 
constantly improved. 

A decade has passed since a fruitful connection be- 
tween wireless communications and statistical mechanics 
was introduced by a seminal work by Tanaka [|[ . On the 
basis of an analogy between the demodulation problem 
of wireless communications and statistical mechanics of 
disordered Ising spin systems, he successfully clarified the 
potential efficiency of a wireless communication scheme 
known as code-division multiple access (CDMA), which 
is employed in the third generation cellular phone sys- 
tems. Later, this analogy was also utilized in developing 
practically feasible and efficient demodulation algorithms 
for CDMA 

In a recent study the connection to statistical me- 
chanics was extended to another wireless-communication 
scheme, namely, orthogonal frequency-division multi- 
plexing (OFDM), which is today employed in the fourth 
generation cellular phones and the latest Wi-Fi systems. 
According to this scheme, the available frequency do- 
main is divided into sub-domains, and data is transmit- 
ted by the subcarriers associated with those sub-domains 
Q . Because of the orthogonality between the subcarri- 
ers, they can be closely placed in the frequency domain, 
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thereby attaining high-rate-data transmission. However, 
in a mobile radio environment, relative movement brings 
about a Doppler spread, which destroys the orthogonal- 
ity between the subcarriers. This destruction leads to 
the occurrence of intercarrier interference (ICI), which 
rapidly deteriorates the bit error rate. An efficient ICT 
cancellation scheme is, therefore, indispensable in use of 
the OFDM scheme in such environments. In Q , a Monte- 
Carlo-based ICI cancellation scheme was developed on 
the basis of mapping an OFDM model to a variant of 
random-field Ising spin systems. Numerical experiments 
indicated that the developed scheme can achieve signifi- 
cantly better performance than existing standard meth- 
ods in terms of bit error rate. However, the cost of com- 
putation, which grows with the square of the number of 
subcarriers, and the technical difficulty in implementing 
electrical circuits prevent the scheme from being practi- 
cally significant. 

The purpose of this study is to develop an approx- 
imate ICI-cancellation scheme for resolving the above- 
mentioned drawbacks. For this purpose, a mean-field 
approximation (MFA) is utilized in conjunction with the 
analogy between OFDM and random field Ising spin sys- 
tems. Naive MFA requires as much computational cost 
as the Monte-Carlo based cancellation scheme. We show 
that utility of translational invariance in the Fourier do- 
main, which is intrinsic in OFDM, makes it possible to 
develop an ICI cancellation scheme whose computational 
cost is proportional to the number of subcarriers. We 
also show that the performance of the developed algo- 
rithm based on MFA is well captured by a discrete map 
of a single variable and that the fact is supported by 
numerical experiments. 

This paper is organized as follows. In section|TTl we ex- 
plain the OFDM model studied in this paper. In section 
IIII1 we propose a MFA-based ICI cancellation scheme. In 
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section HVl we explore the performance of the proposed 
method and its time evolution, and derive the approxi- 
mated expression of the proposed canceller. The deriva- 
tion of the approximated expression is explained in Ap- 
pendix El Finally, section IVl is devoted to the conclusion 
and summary. 



II. MODEL 

When a time sequence of a signal, x = {x t } (t = 
1, • • • ,N), is transmitted in a mobile radio environment, 
the received symbol, y = {yt}, in the multipath channel 
is expressed by 



M 



yt = 



1 S^h ( ^-r 27r M 
y hp exp I v — 1 



N 



x t + rjt- 



(1) 



The time delay of each path is assumed to be zero for 
simplicity, and the channel noise, {?7t}, is independent 
with respect to the time domain. The number of paths 
is M, the amplitude of each path, {h^}, is distributed 
according to the Rayleigh distribution [3, \i 



as 



P(h p ) = h p exp(- 



' 2 h 



(2) 



and {e p } is the Doppler shift, which is assumed to be 
distributed uniformally in the region e p £ [0, e ma x]- With 
this model, the difference between the maximum and 
minimum values of the Doppler shift is significant irre- 
spective of the sign of the Doppler shift. 

By applying the discrete Fourier transformation to 
cq.([T]), the frequency-domain representation of the trans- 
mitted signal, X = {X k }, and the received signal, 
Y = {Yfc}, where k = 0, • • • , N — 1 and N is the number 
of subcarriers, is given by 



N 



(3) 



i=i 



where n k is the discrete Fourier transform of the channel 
noise. We assume that n k is characterized as an additive 
white Gaussian noise (AWGN) of mean zero and vari- 
ance a 2 . The component of N x N matrix W, called a 
frequency-domain matrix, is given by 



Wi 



kl 
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hp sin(7r(Z -k + ep ))e^ l{ - 1 'TT^^- k+t ^ 



MNsm(n(l -k + e p )/N) 



(4) 



where Wki represents the intensity of the interference 
from subcarrier I to k. When {e p } = 0, the matrix is di- 
agnonal, namely Wm = Sk,l, where S is Kronecker's delta, 
and there is no ICI between any subcarriers. In general, 
the frequency-domain matrix has translation symmetry, 
so the value of each component Wu only depends on the 



difference of the indices, namely k— I. This fact is a result 
of the Fourier representation. 

The Doppler shift, e p , is normalized by the frequency 
separation of subcarriers, Af = 1/N, as e p = /d, p /A/, 
where /d, p is the Doppler frequency at p-th path. The 
parameter e p indicates the influence of the Doppler effect 
on ICI with a given alignment of subcarriers. 

The transmitted bits, channel noise, and received bits 
arc represented as complex numbers. Their real and 
imaginary parts are denoted by X R and X 1 , n R and 
n , and and Y 1 , respectively. Hereafter, they are 
represented as vectors consisting of 2N elements: X = 
[X R ,X I ] T , n = [rc^n 1 ] 7 , and Y = [F i? ,F / ] T , respec- 
tively. T denotes the operation of the matrix transpose. 
The corresponding frequency-domain matrix is redefined 

\W R -W-H r> 

as a 2N x 2N matrix, W = ^.j ^ R , where W n 

and W 1 are the real and imaginary parts of the matrix, 
respectively @. 



III. MEAN-FIELD ICI CANCELLER 

The problem of the OFDM system is to recover the 
original signal by canceling out the inter-carrier interfer- 
ence. The Baycsian framework offers various ICI cancel- 
ing strategies on the basis of the posterior distribution, 



P(X\Y) 



P(Y\X)P(X) 
Ex P(Y\X)P(X) 



(5) 



where the likelihood P(Y\X) is given by the distribution 
of the channel noise At the receiving side, it is as- 
sumed that the channel noise is described by the Gauss 
distribution with mean zero and variance a 2 , 

/ 1 \ N / 2 ( (WX-Y) 2 \ 
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and prior distribution P(X) is the uniform distribution. 
The posterior probability can therefore be expressed as 

111 2N 
P(X\Y) = - exp (g E JiiXiXj - E hiXi) } 

i,j i=l 



-PH(X\J.h) 



(7) 



where cr -2 is identified with the "inverse temperature" 
/3, the interaction matrix and the external field are given 
by J = W T W and h = Y T W, respectively, and Z 
corresponds to the partition function. The interaction 
Ji Jal 



matrix, J 



J 2 



Jl 



also has translational invariance; 



the (k, ^-components of TV x TV sub-matrices J\ and J 2 
only depend on k — I. 

The Hamiltonian of the OFDM system, U(X\J,h), 
defined in cq.©, can be regarded as that for a random- 
field spin model. Unlike typical random field models, the 
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random field of the OFDM model is determined by the 
ICI between transmitted bits and the channel's proper- 
ties. In particular, when the number of the paths, M, is 
equal to 1, the value of the off-diagonal element is much 
smaller than that of the diagonal element; hence, the 
system can be regarded as a single-body problem with 
random fields. 

The maximum a posterior probability (MAP) strategy 



update rule is then given as 



X = argmaxP(X|r), 



(8) 



where arg max is the argument giving the maximum value 
of the function, is guaranteed to minimize the block-wise 
error probability. This strategy corresponds to the search 
of the ground state of Hamiltonian H(X\J, h). However, 
the numerical cost of exactly obtaining the maximizer X 
increases exponentially with increasing number of sub- 
carriers. Zero temperature (/3 — > oo) synchronous dy- 
namics of N symbols based on the mean-field approxima- 
tion (MFA) is a practically feasible approximate scheme 
for finding the MAP solution of eq. ((8]) [10| . In the case of 
the quadrature-phase-shift-keying (QPSK) modulation, 
where the components of X takes one of the two values, 
±1, the scheme is expressed as 



X 



1 = sgn(/i fe - E>i*i), 

l^k 



(9) 



where sgn(u) denotes the sign of u and Xt is the ten- 
tative decision of the fc-th symbol after t iterations. 
It is assumed that the configuration {A^;} is invariant 
at {X£} after sufficient updates, and the fixed config- 
uration is regarded as the final decision of the trans- 
mitted bits. Synchronous update schemes similar to 
cq. @ have been introduced for evaluating the minimum 
mean square error (MMSE) estimator of Gaussian priors 
P(X) cx exp (— \X\ 2 / (2(j\)) , in which the transmitted 

bits are estimated as Xmmse = ((c 2 / 'o 2 x )I + J)~\ h m 
the current case, where I is the identity matrix [lll - ll3l |. 

The computational cost of cq.© increases as 0(N 2 ), 
and it may reduce the practical feasibility of eq.©. To 
reduce this cost increase, the proposed algorithm utilizes 
the fact that the absolute value of Ji fc; and Jiu decreases 
as \k — 1\ increases, which indicates that the major contri- 
butions to ICI (to which each subcarrier is subject) are 
concentrated on the subcarrier's neighborhood. Based 
on this observation, the strategy proposed here considers 
only a part of the ICI among subcarriers in the frequency 
domain at each stage of the cancellation. A similar strat- 
egy was also proposed for a Gaussian MMSE estimator 

Let us define the set of indices of subcarriers that are 
considered to be contributed on ICI of subcarrier k as 
d k (uj) = {x,x + N\x = mod(fc + N - 1 ± y, N) + 1, y e 
N, 1 < y < w}, where mod(fc, N) means the remainder of 
k/N as an integer. In general, dk(u) is composed of 4w 
elements; for instance, <9i(l) = {2,N,N + 2,2N}. The 



sgn(h h 



E J ^ x t)- 



(10) 



The algorithm (J9)) corresponds to the case that oj = [N J 
2J in cq. (fTU|) . where [x\ denotes the largest integer not 
greater than real number x. The numerical cost per iter- 
ation of eq. ([T0f is O(N) as long as parameter u is O(l). 
When u < N/2, fixed point {X*} does not correspond 
to the maximizer of the posterior probability, but it is 
expected to be a good approximation of the maximizer. 



IV. RESULTS 



Performance of ICI canceller 



We observe bit error rate (BER), which is defined by 



2N 

BER =2iv"E<W 



(11) 



where 7TT and (• • • ) represent the average over the 
frequency-domain matrix W and over the channel noise 
and transmitted symbol, respectively. The BER perfor- 
mance of the decoder is bounded from below by that for a 
single bit transmitted through the AWGN channel, which 
is given by 



BER' 



OPT 



1 



-erfc(VEb/N0), 



(12) 



because there is no ICI involved when only a single bit is 
transmitted. BER of the ICI canceller closes to BER OPT 
as the elimination of ICI becomes successful. In the cur- 
rent model, the signal to noise ratio (SNR) is given by, 



SNR 



1 

N 



N 

E 



Wi, 



(13) 



In the QPSK modulation case, the number of bits per 
symbol is two, so Eb/NO corresponds to SNR/2 Q. 

We check the performance of the matched filter, whose 
mathematical manipulation corresponds to multiply the 
received symbol by the Hermit conjugate matrix of W . 
The matched filter minimizes the power of the channel 
noise, but it cannot reduce the error due to the intercar- 
rier interference. Therefore, as a canceller of intcrcarrier 
interference, the performance of the proposed algorithm 
should be better than that of the matched filter. In the 
proposed algorithm, uj = corresponds to the demodu- 
lation by the matched filter and as uj increases, the per- 
formance is expected to improve. 

Eb /NO dependence of BER is shown in FigQ] for num- 
bers of the path, (a) M = 3 and (b) M = 15, respec- 
tively. The number of the subcarriers is N = 32, and 
the maximum value of the Doppler shift is e max = 0.5. 
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FIG. 1. (color online) Dependence of BER on Eb/NO at (a) 
TV = 32, M = 3, e max = 0.5 and (b) N = 32, M = 15, e m ax = 
0.5. Results for ui = 1 and uj — 6 are shown and "Full" cor- 
responds to uj = 16. The dashed line and dashed-dotted line 
represent the optimal and matched filter cases, respectively. 
GA means the BER performance under the Gaussian approx- 



imation. 
W. 



The data points are averaged over 10 samples of 



BER OPT and BER for the matched filter are shown by 
the dashed line and dashed-dotted line, respectively. The 
BER performance of the canceller given by eci. tfTU)) with 
u> = 1, which is the simplest case, is better than that of 
the matched filter, and it improves as the number of in- 
teractions increases. When all interactions between the 
subcarriers are taken into account, the BER performance 
of the canceller given by eq. (O almost coincides with the 
optimal performance BER OPT . 

The w-dependence of BER performance at M = 
3, e max = 0.5, and Eb/N0= 8.69 is shown in FigfJ for 
N = 64, 128, and 256. This graph indicates that BER 
does not depend on the number of subcarriers, N. BER 
rapidly decreases as u> increases from zero, and it grad- 
ually approaches BER OPT as oj further increases. At 
uj > 16, the differences between the BER of the proposed 
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FIG. 2. (color online) Dependence of the BER performance 
on u at M = 3, e max = 0.5, and Eb/N0= 8.69dB. Results for 
numbers of subcarriers N = 64, 128, and 256 are shown, and 
the dashed line shows the BER of the optimal limit. 



method and the optimal limit is less than 10% of the 
value of BER OPT . This result indicates that u> ~ 16 is 
sufficient to practically achieve the BER of the original 
MFA based canceller ([9]) irrespective of N; thus, the re- 
quired numerical cost per bit does not increase as O(N). 
In this model, the unit width of the frequency domain 
is given by 1/iV, then the proposed algorithm will pro- 
vide more effective use of the frequency domain without 
increasing the numerical cost. 



B. ICI dynamics 

To determine the validity of the proposed method as 
a realistic canceller, the time evolution of the BER of 
the proposed algorithm, decoding algorithm at N = 
128, M = 15, and e max = 0.5 is plotted in FigE] The 
initial condition, {A^!}, is obtained by the matched filter, 
so BER at the 0-th step corresponds to that at the fixed 
point for uj = 0. The horizontal axis means the number 
of time steps, and the vertical axis represents the BER 
with respect to the tetative decision at the time step. As 
can be seen in the figure, the BER performances at ui = 3 
and oj — 64 (full) converge to equilibrium values after two 
updates. The time steps required to reach the fixed point 
are quite short and only slightly depend on N, and the 
proposed algorithm is useful for an implementation as an 
ICI canceller. 



C. Macroscopic description 

To analyze the time evolution of the canceller given by 
cq. (|10l) . we attempt to describe the dynamics by using 
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a finite number of macroscopic variables [J, The 
overlap between the transmitted bit and the predicted 
bit at step t under a given realization of W is defined as 



2N 

— Y 

2N ^ 

fc=i 



( X k X k 



(14) 



and the BER at step t is given by 0.5 x (1 — to*). The 
simplest description of the time evolution of eq. (|14[) is 
provided by ignoring all spatial/time correlations among 
the subcarrier symbols, which leads to a discrete map of 
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= J Dz sgn(z + yfi%z), (15) 

^ 2N 

-^J fcfc , E? = 2A(1 - to 4 ) + B + <j 2 C, 



2N 
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where Dz = dze z / 2 /y/2~K, and the coefficients for the 
fixed sample of W are given by 
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fc=l l£d k {w) 
^ 2N 

k=n$a k {w) 

2N 
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2N 



(16) 
(17) 
(18) 



fe=i ( 



The derivation is given in detail in Appendix|XJ These ex- 
pression indicate that the ICI is approximated by Gaus- 
sian noise, and noise variance a 2 is effectively increased 
by B/C due to the insufficiency of the ICI cancellation 
of the MFA canceller (fT0|) . The macroscopic equation 
corresponding to cq.® is obtained by setting ui = N/2. 

BER defined at fixed points given by eq. lfTB"]) . denoted 
by GA (Gaussian approximation), are compared to the 
real BER curve in FigiU where BER = (1 - m t )/2. BER 
at the fixed point are in good accordance with the ex- 
perimental data irrespective of the value of u. The time 
evolution of BER is also well described by eq. (p~5|) . as 
shown in FigOJ in which the initial condition, {m®}, is 
chosen to correspond to the BER of the matched filter. 
These results differ substantially from that of the random 
spreading codes, where time correlation plays a signifi- 
cant role in the macroscopic dynamics Q ■ The difference 
implies that the orthogonality between the subcarriers in 
OFDM reduces the time correlation and enables the ICI 
dynamics to converge within a few steps. 



V. CONCLUSION 

A practically feasible ICI canceller for the OFDM 
model, which can be regarded as a variant of the random 
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FIG. 3. (color online) The time-evolution of BER at N = 
128, M = 15, e m ax = 0.5, and Eb/N0=7.45dB. Results for 
to = 3 and full case (cj = 64) are shown, and GA means 
the corresponding Gaussian approximation. The trajectory 
is averaged over 10 samples of W. 



field spin model, was developed. The cancellation scheme 
was derived by applying the mean-field approximation to 
the maximization of the posterior probability which cor- 
responds to the search of the ground state of the spin 
model. The properties of the frequency-domain matrix, 
i.e., translation symmetry and smallness of off-diagonal 
elements compared to diagonal elements, was focused on, 
and only a part of the ICI among the subcarriers located 
within a distance w in the frequency domain for each bit 
was considered. The numerical cost of the ICI canceller 
is thus controlled by the parameter u>. When oj = |^7V/2J , 
the ICI between all subcarriers is considered, and the ICI 
canceller corresponds to the approximated MAP demod- 
ulater. 

The BER of the proposed algorithm used for QPSK 
modulation is better than the matched filter even if u = 
1, and it practically approaches the optimal limit as u> 
increases further. The BER performance is saturated 
near the optimal limit at around a certain uj ~ O(l) in 
the whole Eb/NO region. Furthermore, the performance 
under a given u> only slightly depends on N. This result 
means that the required numerical cost per iteration to 
achieve a feasible performance level is 0(N). 

The number of time steps required to reach a fixed 
point of ICI canceller is 0{\) only slightly depends on 
the value of N and u>. The total numerical cost to elim- 
inate ICI and demodulate the transmitted bits is there- 
fore O(N). The proposed algorithm will be practical to 
implement by virtue of its low computational cost. 

The fixed point of the ICI canceller and the dynamics 
to reach there are well described by a discrete map of 
a single macroscopic variable under the approximation 
of ICI for each bit as independent Gaussian noise. It 
is considered that the orthogonality between the subcar- 
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riers prevents a time correlation being induced, and is 
a mathematical background of the high accuracy of the 
proposed ICI canceller. 

The proposed algorithm is efficient for the parameter 
region where the Doppler shift causes large intercarrier 
interference. By introducing this algorithm as an ICI 
canceller, the OFDM scheme is useful when subcarriers 
are closely arraigned in the frequency domain and a mo- 
bile object moves at high speed. More efficient use of a 
given frequency domain and enhanced accuracy in satel- 
lite communications are also expected. 



any j. With these quantities, the following expression 
can be obtained: 

m\ = J DzBgnlJa + XiJlfz}, (A5) 

where E' 2 = (zf 2 ) — {z\) 2 , and coefficient Xj can be ig- 
nored because the function is invariant against transla- 
tion z — > — z. The approximated expression for the full 
case can be obtained by setting u> = N/2, so <9,;(w) con- 
tains all bits except i. 



ACKNOWLEDGMENTS 



We would like to thank Ido Kanter for his helpful com- 
ments and discussions. This work was supported by JSPS 
Fellow No. 23-4665 (AS) and KAKENHI No. 22300003 
(YK). 



Appendix A: Derivation of macroscopic dynamics 

According to the definition of {hk}, eq.([T0|) can be 
transformed as follows, 

27V 2N 

X* +1 = sgn[- ]T J ij X} + Y, J n X i + H W ^- 

jedi(u>) j=i n=i 

(Al) 

By separating the summation of the first and second term 
into three part; j <G <9i(w), j ^di(uj) and j = i, the time 
evolution of m* can be written as 

m* +1 = (sgn[ J u + Xi J ij( X j- Xt j) 

+ Xi ]T JijXj + Xi^Wnnp]) 

= (sgn[J vl +X lZ l}). (A2) 

The right-hand side of eq. (|A2p depends on the random- 
ness; the transmitted symbol and the channel noise, 
through z\ . The average over the randomness can there- 
fore be replaced by the average over z\ according to an 
appropriate distribution. The distribution of z\ is ap- 
proximated by a Gauss distribution. The first and second 
moments of z\ are given by 

(4) = 0, (A3) 

2JV 

(zf)=2 £ 4(1 -m)) +2 £ 4+<7 2 £w£, 

(A4) 

where it is assumed that Xj = Xj with probability (1 + 
m* )/2, and Xj = -X* with probability (1 - m*-)/2 for 
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